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Abstract. We present in this paper a study on the peculiar supergiant HD 172481. Its spectral type (F2Ia), high galactic 
latitude (& = — 10? 37), circumstellar dust, high radial velocity and moderate metal deficiency ([Fe/H] = — 0.55) confirm the 
post-AGB character of this object. A detailed chemical analysis shows slight but real s-process overabundances, however no 
CNO-enhancement was detected. Furthermore, the spectral energy distribution and the TiO bands in the red part of the spectrum 
reveal a red luminous companion. The luminosity ratio of the hot F type component and this cool M type companion Lf/Lm 
is derived for a reddening of E(B-V) = 0.44 (Lf/Lm ~ 1.8) and indicates that the companion must also be strongly evolved 
and probably evolving along the AGB. Neither our photometric data-set, nor our radial velocity monitoring show evidence 
for orbital variability which may indicate that the period is too large for direct binary interaction. Most interestingly, a strong 
lithium resonance line is detected, which yields an abundance of loge(Li) =3.6. Several explanations for this large lithium 
content are explored. 
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1. Introduction 

Although the research on the late stages of stellar evolution 
gained enormous momentum thanks to the IRAS survey back 
in the beginning of the eighties, a detailed understanding of the 
final stellar evolution in general and of the post-AGB phase in 
particular, is still lacking. Since the evolutionary time for an 
object to cross the HR diagram from the tip of the AGB to the 
planetary nebula phase is short (100-10000 years depending 
on the core-mass; e.g. Blocker 1995) the objects are rare and 
the significance of results on individual objects is difficult to 
extrapolate to the whole AGB stellar evolution theory. 

Moreover, in recent years it became clear that binarity can 
influence the final stages in many ways. This complicates the 
general observational picture since results on disguised bina- 
ries may lead to misinterpretations when confronted with sin- 
gle star evolution theory. Stars like the Red-Rectangle, which 
were longtime considered to be proto-typical post-AGB ob- 
jects, were later discovered to be evolved interacting binaries 
(Waelkens et al. 1996). More recently, ISO results on the de- 
tailed infrared analysis of the circumstellar dust spectra re- 
vealed that also the dust evolution can be quite different in 
binaries than in single stars (e.g. Waters et al. 1998, Molster 
et al. 1999). There is growing observational evidence that in 
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evolved binaries with a wide period distribution, part of the cir- 
cumstellar material is trapped somewhere in the system. This 
has a strong impact on the chemical and dynamical evolution of 
the system (e.g. Van Winckel et al. 1995, 1999; Jura & Kahane 
1999). Thanks to the trapped dust, binaries will have a pro- 
longed IR lifetime so any sample of post-AGB objects defined 
on IR measurements will be biased towards binarity (Waters et 
al. 1997; Van Winckel 1999). 

Since the stellar luminosity increases by large factors (a 
typical L* is between 3000-8000 Lq) during stellar evolution, 
most of the systems with one component in a post-AGB evo- 
lutionary phase, are single-lined spectroscopic binaries (SB1). 
Double-lined binaries should be rare since they must consist of 
two, almost equally luminous and thus almost equally evolved 
objects. The initial mass of both stars must have been therefore 
extremely similar. 

One possible such object was found to be AFGL4106 
(Molster et al. 1999, Van Loon et al. 1999) but it was recog- 
nized to be a massive double evolved system consisting of a 
massive F-type post-red-supergiant and an M-type red super- 
giant. The spectral energy distribution (SED hereafter) is com- 
plex (see Fig. 2 and 5 in Molster et al. 1999) with the F-type 
dominating the UV and optical wavelengths while the M-type 
component dominates the near-IR. In the far-IR the huge cool 
dust component is clearly observed. The luminosity ratio of 
Lir/Ljvf=1.8 was found with an estimated total luminosity of 
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Table 1. Basic parameters of HD 172481 (SIMBAD database). 
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Fig. 1. Geneva photometry (upper panel) and radial velocity 
measurements (lower panel) of HD 172481. 



Lp = 1.3 x 10 5 Lq assuming a distance of 3.3 kpc (Molster et 
al. 1999). 

In this paper we discuss in detail another such rare double 
evolved SB2 system: HD 172481 (Table [l]). The binary consists 
of a luminous F-type star, probably in a post-AGB evolutionary 
stage, and an M-type luminous companion, probably evolving 
on the AGB. The paper is organized as follows: in Sect. 2 we 
give an overview of the divers observational data we acquired 
during the last years. In Sect. 3 we list all the observational ev- 
idence for the binary SB2 nature of the system; we construct 
and discuss the spectral energy distribution (SED) and give the 
spectral indications of the cool component. In Sect 4 we focus 
on the photospherical chemical composition of the F-type com- 
ponent and discuss the high Li content. The overall discussion 
on the system is given in Sect. 5 while the summary of the main 
conclusions is given in Sect. 6. 

2. Data 

2. 1. Photometry 

We obtained 54 photometric datapoints of HD 172481 in the 
Geneva photometric system during several runs (1989-1996) 
within the framework of our photometric monitoring program 
of optically bright post-AGB candidate stars (e.g. Bogaert 
1994). The Geneva V magnitudes are displayed Fig. |] The ob- 
ject is found to be clearly variable with a small peak to peak 
amplitude of 0.22 magn. in V, much larger than the estimated 
observational error of 0.005-0.007 magn (Rufener & Nicolet 
1988). The near-IR JHKLM photometry was taken with the 
ESO lm telescope during the night of 15 to 16/8/92. For the 
construction of the SED (Sect. 3.2) we used the mean magni- 
tudes of the Geneva and near-IR data complemented with pho- 
tometry from Humphreys (1975), Olsen (1984) and the IRAS 
PSC via the SIMBAD database. The Geneva photometric mea- 
surements can be obtained from the authors upon request. 

2.2. Radial Velocities 

Our sample of radial velocities was obtained with the 
CORAVEL radial velocity spectrometer installed at the 
Danish 1.5m telescope at La Silla, Chile. Since the spectrum 
is correlated with a hardware mask of Arcturus (K2III), always 
the same lines contribute to the radial velocity determination 
which makes this sample very homogeneous. The Arcturus 
mask can be used for supergiants with spectral type up to F0 
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and the velocities have an accuracy of ~0.5kms _1 (Baranne 
et al. 1979). Some recent velocity measurements obtained with 
the CORALIE spectrograph mounted on the new Swiss Euler 
telescope at La Silla and also the velocities from our high res- 
olution spectra (Table || for details) were added to this sample. 
All velocity data are plotted in Fig. [j] 

2.3. Spectra 

High resolution, high signal-to-noise spectra were taken with 
the Utrecht Echelle Spectrograph (UES) mounted on the 4.2m 
William Herschel Telescope (WHT) located on La Palma, 
Spain, with the EMMI Spectrograph mounted on the 3.5m 
New Technology Telescope (NTT) located in La Silla, Chile, 
and with the Fiber-fed Extended Range Optical Spectrograph 
(FEROS) installed at the ESO 1 .52-m telescope also at La Silla. 
A log of the observations can be found in Table ^. 

The WHT and NTT spectra were standard reduced using 
the specific echelle package of the MIDAS data analysis system 
and included order tracing, bias correction, cosmic hits identifi- 
cation and cleaning, scattered light subtraction and flat fielding. 
The orders were extracted by using a simple sum of the pixels 
and the wavelength calibration was performed using a ThAr 
spectrum. During the NTT measurement of 1998 the internal 
focus of the instrument camera jumped with a small amount 
of encoder steps, so the full resolution was not obtained. For 
the FEROS spectrum, we used the pipeline reduction which 
includes all steps of a standard echelle reduction. All spectral 
orders were normalized by fitting the continuum through in- 
teractively defined continuum points. Sample spectra can be 
found in Fig. £§,§, § and 

3. HD 172481 as SB2 

3. 1. TiO-bands 

The optical high-resolution spectra confirm the F-spectral type 
given in SIMBAD which is illustrated in Fig. ||. However, in the 
1995 UES spectra, we realized a decreasing S/N with increas- 
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Table 2. Log of the observations (spectra), together with an approximate value for the signal-to-noise of the different spectra, 
the heliocentric radial velocities Whciio and the number of spectral lines used in the velocity determination. The accuracy of 
these velocities is about 1.5 kms , except for the velocities derived from the EMMI spectra (^2.5kms _1 ). The s/n ratio was 
determined by computing the standard deviation in a stellar continuum region. 
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Fig. 2. Comparison of the Hp region of HD 172481 
(WHT+UES) and the massive F supergiant HR 1865 confirm- 
ing the F-type spectrum of HD 172481. 
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ing wavelength, and we attributed this to the increasing contri- 
bution of a cool companion. Our assumption was confirmed by 
the detection of strong red TiO band heads in our 1997 EMMI 
spectra, and later in the 2000 FEROS spectra (Fig. [}]). This is a 
clear indication that a red companion is contributing to the total 
spectrum of HD 172481. Moreover, the strength of these bands 
is clearly variable as evidenced by our 1998 EMMI spectra (see 
also Fig. ||), pointing to a high amplitude variation of the com- 
panion. In the latter spectrum, only weak TiO band heads can 
be seen, but a clear CI multiplet could be identified. Since the 
TiO bands are clearly visible, the cool component must be of 
spectral type M but a more precise spectral classification is dif- 
ficult due to the composite nature of the spectrum. 

3.2. SED 

In Fig. f| the spectral energy distribution (SED) of HD 172481 
is compiled from the photometric data. We used model atmo- 
spheres of Kurucz (Kurucz 1993) to fit the data. We adopted 
the atmospheric parameters of our abundance analysis (Sect. 
0; T c ff = 7250K, log g= 1.5). By calculating synthetic Geneva 
colour indices using the passbands from Rufener & Nicolet 
(1988) and the above Kurucz model, we obtained a total red- 



Table 3. The contribution in percent of the F-type component 
to the total composite spectrum as deduced from the SED fit. 
This information is therefore only indicative and does not take 
into account the variability of the M-type component (see Fig. 
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dening of E(B-V) = 0.44, which is somewhat higher than the 
value published by Bersier (1996): E(B-V) = 0.271. Using a 
single object, the fit of only the F-type model atmosphere 
clearly did not satisfy from about A ~ 10 3 nm redwards 
and the spectral characteristics were unlike dust emission. To 
fit all the data we invoked a second, cool component, with 
T o ff=3500K but the model parameters of the cool compo- 
nent are only derived by a few photometric points and are 
therefore highly uncertain. Since, however, the F component 
is clearly a luminous object, also the M component must be lu- 
minous. The luminosity ratio of both components is found to 



be Li 



F— type /La/ — type ' 



1.8 but this is rather dependent on the 



adopted reddening: a lower reddening of E(B-V) = 0.2 would 
decrease this ratio to 1 . The integration of the dereddened pho- 
tometry gives a luminosity of the F-type component of 600 d 2 
Lq with d in kpc. 

In Table [}] we list the contribution of the F-star to the to- 
tal composite spectrum for some reference wavelengths. From 
about 975 nm redwards, the cool component is more luminous 
than the F-component. 

In addition to the two stellar components, the IRAS fluxes 
(f 12 = 5.41 Jy, / 25 - 5.22 Jy, / 60 - 0.59 Jy, f wo < 
1.85 Jy), which are also plotted on Fig. [ij point to a cool in- 
frared excess which we interpret as caused by a circumstellar 
dust shell or disk. A small excess seems to be already present 
in the L (3.8 /im) and M (4.7 /im) photometric bands and might 
indicate a present day mass-loss of the cool component. Clearly 
more infrared data is needed to characterize the properties of 
the M-type component. 
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Fig. 3. The EMMI and FEROS spectra in the 700-715 nm TiO band head region. The spectra are offset by 0.3 units. The TiO 
bands in the red spectrum of HD 172481 are a clear signature of its cool companion. Moreover, these bands vary in time as 
evidenced by our 1998 EMMI spectrum. Band head wavelengths are from Valenti et al. (1998). 
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Fig. 4. The dereddened SED of HD 172481 with a E(B-V) 
= 0.44. The dotted line represents the Kurucz model with 
T c ff=7250K and log(g)=1.5; the dashed line the model with 
T o ff=3500K and log(g)=0.5; the full line the combined model. 
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We can conclude that the SED consists of 3 components: a 
luminous F component dominating the UV and visible part of 



the spectrum, a luminous M component dominating the near-IR 
and a cool dust envelope dominating the IRAS fluxes. 

3.3. Photometry and Radial Velocity 

It is clear from Fig. |that HD 172481 is variable, in magnitude 
as well as in radial velocity. Despite the poor sampling, sev- 
eral techniques (CLEAN Roberts et al. 1987, PDM Stellingwerf 
1978, and SCARGLE Scargle 1982, Home & Baliunas 1986) 
were used to find a period (pulsational or orbital) in these data, 
but none of them yielded a convincing stable period, in the 
sense that the folded light curve never satisfied. Nevertheless, 
by fitting details of the light curve as has been done in Fig. 
g, we can conclude that HD 172481 is a semi -regular pulsator 
varying with a timescale of about ^100 days. 

One can argue that the variability in the photometry is in- 
duced by a variable contribution of the M-type companion in 
the Geneva V filter. However, this seems very unlikely be- 
cause the latter filter is centered around 549 nm where the 
contribution of the cool companion is only ~10% (Table |]). 
Moreover, this filter correlates well with the U photometric fil- 
ter (A c = 346 nm) for which the contribution is negligable. Even 
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in the G filter ( A c = 58 1 nm), no clear trace of a variable contri- 
bution could be found. In addition, one would expect a larger 
amplitude if the pulsation is (also) due to the cool AGB star: an 
amplitude of 0.2 mag. is rather usual for post- AGB stars. 

In Fig. |^ we present our high resolution spectra around the 
Ha line. The line exhibits a shell type profile, an emission pro- 
file that is found in almost all other optically bright post- AGB 
candidate stars (e.g. Waters et al. 1993; Oudmaijer & Bakker 
1994; Van Winckel et al. 1996). It is still unclear which mech- 
anism is responsible for the variable Ha emission in the spec- 
tra of post- AGB stars. A recent photometric and spectroscopic 
study about the nature of the pulsation of the bright post- AGB 
star HD 56126 (Barthes et al. 2000) suggests that the irregular 
pulsation probably originates from the interaction between the 
photospheric dynamics and the dynamics of the higher atmo- 
spheric layers. Shock waves are generated when the receding 
photospheric layer collides with the upcoming one and proba- 
bly play a major role both in driving the pulsation and limiting 
its amplitude. 

Although a peak-to-peak amplitude of 17kms _1 is ob- 
served, we could not discover a long-term trend in the radial 
velocity data which would point to resolved binary motion. 
This could be a projection effect, if we view the system on 
a unfavourable inclination, but is likely due to the long orbital 
period and hence very small velocity amplitude. Assuming typ- 
ical masses for the post-AGB and AGB objects (F-type = 0.6 
Mq and M-type =1.0 Mq), an inclination of 60° and a veloc- 
ity amplitude of 5 km s ~ 1 , we obtain a period of 5 1 years . Note 
that the CORAVEL mask is defined in the blue (A < 520 nm) 
so the contribution of the M-component in negligable. 

We can conclude that both the light curve and the radial ve- 
locity data give evidence for the complex and not well under- 
stood variability of the atmospheric layers seen in most post- 
AGB stars. This is further strengthened by the variable Ha 
emission. Evidence for binary induced photometric variability 
or orbital motion is not observed, which is probably indicative 
of a long orbital period. The radial velocity sample indicates 
too long an orbital period to expect direct interaction during 
stellar evolution. 

4. Chemical Analysis 

We performed a detailed abundance analysis on our UES spec- 
tra to gain insight into the evolutionary status of the F-type 
component. We confined our analysis to these two spectra for 
several reasons: due to the increasing contribution of the cool 
companion at longer wavelengths, we decided not to use the 
EMMI spectra since they concentrate on the red region of the 
optical spectrum. Since many useful lines are situated in the 
blue and green region a significant chemical study can still be 
performed on this limited spectral domain. Moreover, the UES 
spectra display a very good signal-to-noise. 

Apart from the resolution difference, the quality of the 
FEROS spectrum is very similar to the UES spectra (Fig. 0), 
but a comparison of the red region reveals that the contribu- 
tion of the cool companion was larger during the FEROS run. 
Therefore, we preferred to use solely our UES spectra for the 
chemical analysis. 



5 

Fig. 5. Detail of the light curve of HD 172481. The harmonic 
fit with a period of P=91.7d is only considered to indicate the 
typical timescale of the semiregular variability of HD 172481. 
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4. 1. Line data and Atmospheric Parameters 

In order not to compromise the results of a chemical analy- 
sis, one has to take care that only lines with well determined 
oscillator strengths (log ^/-values) are selected. A list of such 
lines useful for the chemical analysis of A and F type stars has 
been collected at our institute during the past few years and is 
still regularly updated. We described this list in detail in Van 
Winckel & Reyniers (2000). The list we used for HD 17248 1 is 
essentially the same since only a few lines were added. 

We used the CDROM-grid of LTE model atmospheres by 
Kurucz (1993) in combination with his abundance calculation 
programme WIDTH9. The parameters of an atmospheric model 
are determined by forcing the computed abundance of Fe to be 
independent of excitation potential (determination of T c fr), ion- 
isation stage (gravity) and reduced equivalent width (microtur- 
bulent velocity). This method yielded T ff = 7250 K, log g = 1 .5 
and £ t = 4kms~ 1 . 48 Fel-lines and 14 Fell-lines were used 
to determine the appropriate atmospheric parameters. With this 
model, we also obtain a very good abundance consistency for 
all other chemical elements for which lines of two ionisation 
stages are available (Ca, Cr, Mn and Ni) confirming an accu- 
rate set of atmospheric parameters. 

The equivalent widths of 209 lines (28 different ions) used 
in the analysis were measured, using direct integration for un- 
blended lines and multiple Gaussian fitting for blended lines. A 
rather surprising discovery in the spectra we took is the pres- 
ence of the strong Li resonance line at 670.8 nm (see Fig. ^| and 
. Before discussing in detail the overall results of the chemi- 
cal analysis we focus on the determination of the Li abundance. 

4.2. Li-abundance 

The lithium abundance is derived by spectrum synthesis of the 
670.8 nm region. Spectrum synthesis is used to account for the 
doublet structure of the Li resonance line. Line parameters are 
taken from Kurucz (1993) and are nearly the same as in Cunha 
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Fig. 7. Atomic lines of s-process elements and Li I and bands of TiO in one spectrum. On the left panel, lines of s-process 
elements can be seen, with the remarkable detection of a Hf II line on our UES spectra. On the right panel, the 671 nm region is 
plotted, with the Li I doublet. Clearly, the contribution of the cool companion in this region is larger for our FEROS spectrum 
(which is offset by 0.3 units) as evidenced by the detection of the TiO band head at 67 1 .447 nm. 
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et al. (1995). Kurucz's program SYNTHE was used to make the 
syntheses. 

Synthetic spectra have to be convolved with a macroturbu- 
lent broadening function and with the instrumental broaden- 
ing, to match the line widths of the observed spectrum. For 
the instrumental broadening we used a gaussian with a full 
width at half maximum given by the spectral resolution of the 
UES spectrograph: Cinstr = 6kms~ 1 . The macroturbulent ve- 
locity Cmaci is determined in the following way. We selected 
17 clearly unblended lines of different width (< 150 mA) and 
different species. For each line, we calculated the abundance 
from the equivalent width using the WlDTH9-program. After 
that, we synthesized with the Kurucz spectral synthesis pro- 
gramme SYNTHE each line using the obtained abundance. The 
total broadening velocity Ctot (Ctot = CLtr+CLcr) is the onl y 
free parameter in this synthesis. We found a slight dependence 
of Ctot on the equivalent width of the line: 

Ctot = 0.09* W x + 9.18 

with a classical correlation coefficient of 0.93. We adopted for 
Ctot a value of 14kms _1 (hence Cmacr — 12.6kms _1 ) for the 



Table 4. Sensitivity of the Li-abundance to the stellar param- 
eters and an upper limit for the non-LTE correction. The stan- 
dard parameters of HD 172481 are T cff =7250K, log g = 1.5, 
Ct = 4 km s -1 and C = 8kms _1 ; the standard Li-abundance 
loge(Li) = 3.57. 
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* log e(Li)LTE— log e(Li) non _LTE, estimated from Carlsson et al. 
(1994). 

Lil-line, yielding a Li-abundance of loge(Li) — 3.57 using a 
least-squares fit (Fig. |8j). 

In Table [| we investigate the sensitivity of the derived Li- 
abundance to the uncertainties in the stellar parameters. It is 
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Fig. 6. Normalized spectra of HD 172481 around Ha. Each 
spectrum is offset by 1 unit. 
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clear that the temperature is by far the most delicate parame- 
ter, though a variation of 250K (as in Table is rather large. 
We estimate the uncertainty of the lithium abundance around 
0.20 dex. To our knowledge, detailed non-LTE calculations of 
the Li I resonance line for this temperature-gravity domain are 
not available in literature, but an extrapolation of the results 
published by Carlsson et al. (1994) indicates that a non-LTE 
correction would not exceed 0.2 dex (i.e. < loge(Li)LTE — 
log e(Li) non _ LT E < 0.2). 

Adopting loge(Li) — 3.57, we expect for the sub- 
ordinate Lil-line at 610.36 nm an equivalent width of 
W\ = 5 mA. Unfortunately, this line is blended with a Fell- 
line (A610.350nm). Nevertheless, an upper limit for the Li- 
abundance can be derived from this line: log e (Li) < 3.7. 
The Lil-line at 8 12.645 nm has a calculated equivalent width 
of 0.6 mA and is obviously not detectable on our spectra. 

4.3. Overall results of the chemical analysis 

The derived abundances of all other elements are summarized 
in Table ||; a graphical representation can be found in Fig. ^. 
Besides the lithium line at 670.8 nm, also the resonance line of 
europium (A 412.973 nm) was treated with spectrum synthesis, 



Fig. 8. The Lil resonance line in the UES-spectrum of 
HD 172481. The points are the observed spectrum, the lines 
synthetic spectra with loge(Li) = 3.47, 3.57 and 3.67 resp. A 
least-squares fit resulted in log e(Li) = 3.57 (full line). 
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in order to account for the strong hyperfine splitting of this line. 
The hyperfine data were taken from Biehl (1976). 

The carbon (C) and oxygen (O) abundances are well es- 
tablished, using the multiplets at 615.6 nm and 645.4 nm for 
oxygen. Carbon is clearly not enhanced, and the slight over- 
abundance of oxygen ([O/Fe]=+0.2) is normal in this metal- 
licity range (e.g. Boesgaard et al. 1999). A nitrogen (N) abun- 
dance was hard to obtain, because most useful N lines are sit- 
uated in the red part of the spectrum, which we did not use in 
our analysis. The N abundance is derived from very small lines 
(equivalent widths 11, 5 and 7 mA). Consequently, one has to 
interpret the tabulated abundance (log e(N) = 7.61) more as an 
upper limit rather than as a fixed abundance. We can conclude 
that also N is not enhanced. 

No sodium (Na) enhancement is measured. The a-elements 
Mg, Si, S, Ca and Ti are all enhanced with respect to iron, 
with a simple mean [a/Fe]=0.36. Such a value is expected 
in this metallicity range (e.g. McWilliam 1997). The titanium 
(Ti) abundance is quite high, but it is deduced from rather 
strong lines (7 out of the 13 lines have W\> 100mA) and 
therefore very sensitive to the micro turbulent velocity: an in- 
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crease of lkms -1 in £ t decreases the Ti abundance by 0.1 
dex. Scandium (Sc) is also quite high, but this element behaves 
like an even-Z a-element (Nissen et al. 2000) and it is also 
known to be slightly hyperfine structure sensitive. The high 
vanadium (V) abundance is less easy to explain, because nor- 
mally it should follow the iron deficiency (Chen et al. 2000). 
The hyperfine structure sensitivity alone cannot explain this 
overabundance. We searched for other V II lines from the list of 
Biemont et al. (1989), but, using the extensive lists of Kurucz 
(1993) and the VALD-2 database (Kupkaet al. 1999), all turned 
out to be blended. The iron peak elements (Cr, Ni, Zn) do fol- 
low the iron deficiency, only manganese (Mn) is slightly under- 
abundant, which is also a Galactic evolution effect (Nissen et 
al. 2000). 

Moderate but clear s-process overabundances are found for 
HD 172481. The strontium (Sr) abundance is very difficult to 
determine accurately since the resonance lines at A421.5nm 
and A407.7nm are heavily saturated and weaker optical lines 
are not present. The other two elements of the Sr-peak around 
the magic neutron number 50, yttrium (Y) and zirconium (Zr), 
have more usable lines and turn out to be clearly enhanced by 
0.5 dex. For the heavy s-process elements around magic neu- 
tron number 82, the situation is less clear, due to the smaller 
amount of lines. Nevertheless, we feel confident that also the hs 
elements lanthanum (La) and cerium (Ce) are enhanced. Just 
like Sr, an accurate barium (Ba) abundance determination is 
difficult due to the lack of weak lines. For neodymium (Nd), 
we deduced an upper limit of loge(Nd) < 0.9 ([Nd/Fe]<0.0) 
from the non-detection of the line at 529.3 16 nm. We did the 
same for samarium (Sm) and obtained loge(Sm) < 0.55 
([Sm/Fe]<0.1) from the non detection of the 453.794 nm line. 
The r-process element europium (Eu) is slightly overabundant 
relative to iron, but this behaviour is expected for this metallic - 
ity (e.g. Woolf et al. 1995). So Nd, Sm and Eu seem not to be 
intrinsically enhanced. 

Surprisingly, we also detected a small line at A409.316nm 
which we identified as Hf II (Fig. |). The same line is detected 
in the heavily s-process enriched post-AGB stars (Van Winckel 
& Reyniers 2000). An abundance calculation with W\=\2 mA 
yields a very high [Hf/Fe]=+0.6! This line is hardly detectable 
on our FEROS spectrum (also Fig. 0) and the abundance is ex- 
treme but we could not find possible blending lines. Other pos- 
sibly detectable Hf lines are all blended. The significance of the 
high hafnium abundance should be further investigated. 

5. Discussion 

5. 1. Post-AGB nature 

In this section we list the arguments for the proposed post-AGB 
character of our object. 

1. The high galactic latitude (Tab. J} b= -10? 37), 

2. the large radial velocities (Fig. nb and 

3. the moderate metal deficiency ([Fe/H]=— 0.55) 

point to a population II membership. Furthermore, 

4. the SED (Fig. shows the presence of dust, as a result of 
previous mass loss probably during the AGB phase. 



Table 5. Chemical analysis of HD 172481. For every ion we 
list the number of lines used N, the mean equivalent width 
W\, the absolute abundance loge (i.e. relative to H: loge = 
logX/H + 12), the abundance ratio relative to iron [el/Fe] and 
the internal scatter cr, if more than one line is used. For the solar 
iron abundance we used the meteoric iron abundance of 7.51. 
For the solar C, N and O abundances we adopted resp. 8.57, 
7.99 and 8.86 (C: Biemont et al. 1993, N: Hibbert et al. 1991, 
O: Biemont et al. 1991); the other solar abundances are taken 
from Grevesse (1989). ss stands for spectrum synthesis. 

HD 172481 

Tefi = 7250 K 
log g = 1.5 (cgs) 
C i = 4.0kms- 1 
[Fe/H] = -0.55 



ion 


JN 


H A 


loge 


[el/rej 


a 


Li I 


1 


ss 


3.57 




0.20 


CI 


11 


36 


8.01 


-0.01 


0.11 


NI 


3 


8 


7.61 


+0.17 


0.03 


OI 


5 


26 


8.47 


+0.16 


0.14 


Nal 


2 


32 


5.76 


-0.02 


0.12 


Mgl 


1 


37 


7.30 


+0.27 




Si I 


5 


16 


7.41 


+0.41 


0.12 


SI 


4 


32 


7.06 


+0.40 


0.09 


Cal 


15 


55 


5.99 


+0.18 


0.14 


Call 


2 


30 


5.92 


+0.11 


0.02 


Sell 


4 


85 


2.97 


+0.42 


0.03 


Tin 


13 


104 


4.97 


+0.53 


0.15 


VII 


1 


61 


3.85 


+0.40 




QT 


8 


21 


5.17 


+0.05 


0.11 


CrII 


18 


61 


5.17 


+0.05 


0.13 


Mnl 


3 


21 


4.62 


-0.22 


0.01 


Mnll 


1 


13 


4.66 


-0.18 




Fel 


48 


62 


6.92 


-0.04 


0.16 


Fell 


14 


68 


6.96 




0.15 


Nil 


18 


27 


5.77 


+0.07 


0.16 


Nill 


3 


47 


5.78 


+0.08 


0.19 


Znl 


2 


32 


4.17 


+0.12 


0.03 


YII 


6 


42 


2.17 


+0.48 


0.06 


Zrll 


9 


50 


2.55 


+0.50 


0.13 


La II 


4 


70 


1.27 


+0.60 


0.04 


Cell 


6 


21 


1.24 


+0.24 


0.10 


Ndll 






<0.9 


<+0.0 




Smll 






<0.55 


<+0.1 




EuII 


1 


ss 


0.19 


+0.23 




Hfll 


1 


12 


0.91 


+0.58 





5. The photometry and the Ha profiles show a variability very 
similar to other post-AGB stars (Sect. 3.3). 

6. Finally, the slight s-process overabundances are probably 
the result of the third dredge-up, when He burning products 
are brought to the surface. 



The last argument is, however, not straightforward. If we expect 
the signature of a 3rd dredge-up, also the CNO elements should 
be enhanced, which is clearly not detected: the total CNO abun- 
dance is equal to the initial value. 
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Fig. 9. The abundances of HD 172481 relative to iron [el/Fe]. 
Errorbars are plotted if the number of lines is more than 3. The 
horizontal line represents [el/Fe]=0. Details of the adopted so- 
lar abundances can be found in the caption of Table pi 
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Binary post-AGB stars are often depleted, with a dust disk 
as an ideal site for the gas-dust separation (see Van Winckel 
1999 for a review). Depletion alters the abundances profoundly. 
An ideal element to trace for possible depletion is zinc (Zn) 
with its low condensation temperature. A high [Zn/Fe] ratio 
points to efficient depletion. For HD 172481, [Zn/Fe]=+0.12, 
suggesting that no depletion process has taken place. 

5.2. Li-rich evolved stars 

Plez et al. (1993) published an abundance analysis of seven lu- 
minous AGB stars in the Small Magellanic Cloud. These mod- 
erately metal deficient stars ([Fe/H]=— 0.5) are Li rich, C poor 
and show a low 12 C/ 13 C ratio. The chemical pattern observed 
in these stars is explained by envelope burning, which is the 
activation of the H burning CN cycle at the base of the con- 
vective envelope (called "hot bottom burning" HBB). The Li 
abundance in these stars is typically loge(Li) = 3.0, which 
is much larger than expected from standard AGB evolution 
where Li has been destroyed and/or diluted due to convec- 
tion. The mechanism thought to be responsible for this Li pro- 
duction is the 7 Be-transport mechanism (Cameron & Fowler 
1971) in which 3 He is converted to 7 Li via the reaction chain 
3 He( 4 He,7) 7 Be(e, is) 7 Li. Theoretical models suggest that this 
mechanism can only be efficient in massive AGB stars (4-7 
M Q ; e.g. Sackmann & Boothroyd 1992, Mazzitelli et al. 1999). 
Also massive galactic AGB stars with signs of envelope burn- 
ing were found (Garcfa-Lario et al. 1999), but not yet exten- 
sively chemically analysed. 

On the other hand, Li is found in many evolved stars in 
the AGB or post-AGB evolutionary phase for which there is 
strong observational evidence that many of these stars are low- 
mass objects with M<2— 3 M Q : galactic low-mass C stars have 
been found with Li excesses of log e (Li) > 1 (Abia et al. 
1993, Abia & Isern 1996, 1997, 2000); unexplained high Li 
abundances were also found in some metal poor population 
II post-AGB stars, which showed strong evidence of an ef- 



ficient 3rd dredge-up: IRAS 22272+5435 (Zacs et al. 1995), 
IRAS 05341+0852 (Reddy et al. 1997), IRAS Z02229+6208 
and IRAS 07430+1 115 (Reddy et al. 1999). It has been sug- 
gested by Reddy et al. (1999) that these stars could be the suc- 
cessors of the Li-rich carbon stars studied by Abia (see ref- 
erences above). However, a detailed abundance analysis of 12 
J-type Li-enhanced carbon stars (Abia & Isern 2000) showed 
no s-process enhancements in these stars and a solar metallic- 
ity while the 12 C/ 13 C ratio is low. These J-type carbon stars are 
clearly chemically different from the 4 Li-enhanced metal-poor 
post-AGB stars with high 12 C/ 13 C ratios and strong s-process 
overabundances. 

It is generally known that also among red giants ascend- 
ing the 1st giant branch high Li-abundances can be found. 
Wasserburg et al. (1995) developed a non-standard mixing 
mechanism for these RGB stars to explain their low 12 C/ 13 C 
ratio (dubbed "cool bottom processing" CBP) and Sackmann 
& Boothroyd (1999) showed that CBP can account for the 
high lithium abundances observed in the RGB Li-rich K gi- 
ants (de la Reza et al. 1997 and references therein). Moreover, 
they suggested that CBP might also occur in the TP- AGB phase 
(Boothroyd & Sackmann 1999). de La Reza et al. (1997) sug- 
gested that all K giants become Li-rich for a short time. During 
this period the giants are associated with an expanding thin 
circumstellar shell supposedly triggered by an abrupt inter- 
nal mixing mechanism resulting in a surface Li enrichment. 
However, Jasniewicz et al. (1999) did not find any correlation 
between Li abundance and infrared excess in a sample of 29 
late-type RGB giants with infrared excess. There are even Li- 
rich giants with no circumstellar dust. 

Alternative external scenarios to explain the high lithium 
abundances invoke the accretion of a substellar object (objects). 
Such scenarios gained interest by the discovery of several ex- 
tra solar planets and brown dwarf candidates around solar-type 
stars. Planets could be even commonly present around PNe 
(Soker 1999). Engulfment scenarios were developed (e.g. Siess 
& Livio 1999a,b), but cannot account for Li abundances ex- 
ceeding the initial one and therefore fail to explain the very Li 
rich stars with loge(Li) > 2.8 (Siess & Livio 1999b). In other 
scenarios (Denissenkov & Weiss 2000) planet engulfment acts 
as the activating mechanism for the CBP, producing also these 
higher Li abundances. 

We can conclude that the Li-enhancements in the low-mass 
C stars and the 4 metal-poor post-AGB giants emphasize our 
poor knowledge of the dredge-up and mixing processes or pos- 
sible accretion during stellar evolution. 

5.3. Li-enhancement in HD 172481 

It is clear that planet engulfment (alone) cannot explain a Li 
abundance of log e(Li) = 3.57. "Hot bottom burning" (HBB) is 
thought to be responsible for the Li enhancements in moder- 
ately metal-poor luminous AGB stars (estimated initial mass 
~7M Q ) in the Small Magellanic Cloud (Plez et al. 1993). 
This scenario is theoretically thought to occur only in luminous 
AGB stars of initial mass 3-7 Mq (Sackmann & Boothroyd 
1992). The low metallicity and the kinematics of HD 172481 
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suggest an old and hence lower mass nature. Another argument 
against the HBB scenario for HD 172481 is the absence of a 
nitrogen enhancement. During the HBB, C is converted into 
N by the CN cycle. But possibly lithium production may have 
occurred without burning significant amounts of carbon. 

On the other hand, the CBP mechanism may be relevant for 
the Li enhancement in HD 17248 1 . The presence of circumstel- 
lar dust is in favour of de La Reza's scenario, if applicable to 
low mass AGB stars. Unfortunately, we are not able to derive a 
12 C/ 13 C ratio, which could shed some light on the deep mixing 
involved. 

The recent surprising results of the chemical content of J- 
type carbon stars (Abia & Isern 2000) are interesting in many 
respects. The O-rich circumstellar dust which is observed in 
about 5%-10% of these objects (Lloyd Evans 1991) is diffi- 
cult to reconcile with standard AGB evolution and led Little- 
Marenin (1986) and Lloyd Evans (1990, 1991) to conclude that 
the material expelled by the M-type progenitor was trapped 
and stored in a binary system. Recent CO measurements of 
J-type carbon stars confirm this scenario and point to a long- 
lived reservoir of oxygen rich dust (Jura & Kahane 1999). In 
the scenario proposed by Yamamura et al. (2000) to explain the 
ISO-SWS spectra of some J-type carbon stars, these stars are 
wide binaries in which the oxygen rich dust was stored around 
the unseen companion. Note, however, that no direct evidence 
of the binary nature of these stars exists yet. Since no large am- 
plitude radial velocity variations are observed (Barnbaum et al. 
1991), the eventual orbits must be rather wide and direct bi- 
nary interaction by Roche-Lobe overflow is unlikely to have 
happened during evolution on the AGB. Whether the special 
chemical composition of these stars is linked to the proposed 
binarity in wide systems is still an open question. 

For HD 172481 there is direct observational evidence for 
binarity in a probably wide system. Also this object shows 
very peculiar abundances with a high Li abundance, only very 
small overabundances of s-process elements and no strong C 
enhancement. Here, again, the eventual role of the binary na- 
ture of HD 172481 with respect to the remarkable chemical 
composition remains an open question. 

6. Summary 

We presented a study of the peculiar supergiant HD 172481, 
revealing it to be a double lined spectroscopic binary consist- 
ing of an F-type post- AGB star and an M type, probably AGB, 
companion. The post-AGB nature of the F type component is 
deduced from its spectral type; high galactic latitude; circum- 
stellar dust; high radial velocity; moderate metal deficiency and 
slight overabundances of the s-process elements Y, Zr, La and 
Ce. The photometry and the Ha profiles show a variability very 
similar to other post-AGB stars. The spectral energy distribu- 
tion and the TiO bands in the red part of the spectrum reveal bi- 
narity with a luminous M-type companion with a luminosity ra- 
tio of Li?/LM=l-8 assuming a total reddening of E(B-V) = 0.44. 
The lack of a resolved long-term trend in the radial velocity 
data suggests a wide orbit. Most interesting, we found a very 
high lithium abundance of loge(Li) = 3.57. This abundance is 
difficult to reconcile with a "hot bottom burning" scenario be- 



cause of the probably low mass nature of HD 172481. An alter- 
native explanation can be "cool bottom processing" (Sackmann 
& Boothroyd 1999). 
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